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We consider random variables of the form F = fiVi, . . . ,Vn), 
where / is a smooth function and , i £ N, are random variables with 
absolutely continuous law Pi(y) dy. We assume that pi, i = 1, . . . ,n, 
are piecewise differentiable and we develop a differential calculus of 
Malliavin type based on dlnpi. This allows us to establish an in- 
tegration by parts formula E{di4>{F)G) = E(<t>(F)Hi{F,G)), where 
Hi{F,G) is a random variable constructed using the differential op- 
erators acting on F and G. We use this formula in order to give 
numerical algorithms for sensitivity computations in a model driven 
by a Levy process. 

1. Introduction. In recent years, following the pioneering papers [12, 13], 
much work concerning numerical applications of stochastic variational cal- 
culus (Malliavin calculus) has been carried out. This mainly concerns ap- 
plications in mathematical finance: computation of conditional expectations 
(which appear in, e.g., American option pricing) and of sensitivities (the so- 
called Greeks). The models at hand are usually log-normal type diffusions 
and one may then use standard Malliavin calculus. Currently, there is in- 
creasing interest in jump-type diffusions (see, e.g., [7]) and one must then use 
the stochastic variational calculus corresponding to Poisson point processes. 
Such a calculus has already been developed (in [4] and [15]) concerning the 
noise coming from the amplitudes of the jumps and (in [6, 9, 18, 19, 21, 22]) 
concerning times. Recently, Bouleau (see [5]) established the so-called error 
calculus based on the Dirichlet forms language and showed that the ap- 
proaches in [4] and [6] fit into this framework. Another point of view, based 
on chaos decomposition, may be found in [3, 10, 16, 17, 23]. 
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Let us finally mention that the models considered in mathematical finance 
(e.g., Merton's model) may have both a diffusion component (driven by a 
Brownian motion) and a jump part (driven by a compound Poisson process). 
In this case, one may use the standard Malliavin calculus with respect to the 
Brownian motion increments after conditioning in a clear way with respect 
to the Poisson component. This is done in [8, 11, 20]. 

The aim of this paper is to give a concrete application of the Malliavin 
calculus approach to sensitivity computations (Greeks) for pure jump dif- 
fusion models. We give three examples: in the first, we use the Malliavin 
calculus with respect to the jump amplitudes and in the second, we differ- 
entiate with respect to the jump times. In the third, we differentiate with 
respect to both of them. 

The basic tool is an integration by parts formula which is analogous to the 
one in the standard Malliavin calculus on the Wiener space. Here, we give an 
abstract approach which, in particular, permits us to treat in an unified way 
the derivatives with respect to the times and the amplitudes of the jumps 
of Levy processes. More precisely, we consider functionals of a finite number 
of random variables Vi,i = l,...,n. The only assumption is that for each 
i = 1, . . . ,n, the conditional law of Vi (with respect to Vj, j i) is absolutely 
continuous with respect to the Lebesgue measure and the conditional den- 
sity Pi =pi{uj,y) is piecewise differentiable. Using integration by parts, one 
may settle the duality relation which represents the starting point in Malli- 
avin calculus. However, some border terms will appear corresponding to the 
points at which pi is not continuous: for example, if Vi has a uniform condi- 
tional law on [0, 1], the density is Pi{oj., y) = l[o,i] (y) and integration by parts 
produces border terms in and in 1. There is a simple idea which permits 
us to cancel the border terms: we introduce in our calculus some weights 
TTj which are null at the points of singularity of pi — in the previous exam- 
ple, we may take TTi{y) = y"(l — y)°^ with some a G (0, 1). We then obtain a 
standard duality relation between the Malliavin derivative and the Skorohod 
integral and the machinery established in the Malliavin calculus produces 
an integration by parts formula. But there is a difficulty hidden in this pro- 
cess: the differential operators involve the weights vTj and their derivatives. 
In the previous example, i:[{ijj,y) = a(y"~^(l — y)" — y"(l — y)°^~^)- These 
derivatives blow up in the neighborhood of the singularity points and this 
produces some nontrivial integrability problems. We must therefore search 
an equilibrium between the speed of convergence to zero and the speed with 
which the derivatives of the weights blow up in the singularity points. This 
leads to a nondegeneracy condition which involves the weights and their 
derivatives. 

The integration by parts formula is established in Section 2. Since nu- 
merical algorithms involve only functions of a finite number of variables, we 
do not develop here an infinite-dimensional Malliavin calculus, but restrict 
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ourselves to simple functionals. In Section 3, we use the integration by parts 
formula in order to compute the Delta (derivative with respect to the ini- 
tial condition) for European options based on an asset which follows a pure 
jump diffusion equation and in Section 4, we give numerical results. 

2. Malliavin calculus for simple functionals. 

2.1. The frame. We consider a probability space {i},T,P), a sub-ci- 
algebra Q C and a sequence of random variables F^jZ € N. We define 

= ^ V a{Vj,j ^ i). Our aim is to establish an integration by parts for- 
mula for functionals of V^, i S N, which is analogous to the one in Malliavin 
calculus. The cr-algebra Q appears in order to describe all the randomness 
which is not involved in the differential calculus. 

We will work on some set A ^ Q which will be fixed throughout this 
section. We denote by L(oo)(A) the space of the random variables such that 
E(|F|p1a) < oo for all p € N, and L(p^)(y4) will be the space of the random 
variables F for which there exists some 5 > such that E(|Fp+^l^) < oo. 
We assume the following. 

Hypothesis 2.1. € L(oo)(^), i e N. 

For each i € N, we consider some /cj S N and some t/j-measurable random 
variables 

We define 

j=0 

Notice that we may take 

Oj = — oo and bi = oo. 

We will work with functions defined on (aj(a;), which are smooth 
except for the points tj,j = l,...,fcj. We define Ck{Bi) to be the set of 
measurable functions f :Q. x M — M such that for every iv, y ^ f{uj,y) is 
/c-times differentiable on Bi(uj) and for each j = 1, . . . , /cj, the left-hand side 
and the right-hand side limits f{uj,tl{uj)—),f{LO,tl{uj)+) exist and are finite 
[for j = (resp., j = ki + 1), we assume that the right-hand side (resp., the 
left-hand side) limit exists and is finite]. We define 

r.(/) = i:(/(a;,t|M-)-/(^,iiM+)) 

(2.1) 

+ f{u;,bi{u;)-) - f{uj,ai{uj)+). 
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For f,g(z Ci{Bi), the integration by parts formula gives 
(2.2) / fg'{uj,y)dy = T,{fg)- f f'g{uj,y)dy, 

J{ai,bi) J{ai,bi) 

SO Fj represents the contribution of the border terms — or, in other words, 
of the singularities of / or g. 

Let n, /c S N. We denote by Cn^k the class of the G x ;B(M")-measurable 
functions / : $7 x M" ^ M such that /«(/) € Ck{Bi), z = 1, . . . , n, where 

h{f){^, y) ■■= /(w, Vi, . . . , Vi-i,y, Vi+u • • • , K)- 
For a multi-index a = (qi, . . . , a^) G {1, . . . , n}^, we define 

d^f = f 

Moreover, we denote by Cn,k{^) the space of functions / € C„^fe such that 
for every <p < k and every a = (ai, . . . , Op) € {1, ... , n}^, d^f{Vi, . . . , y„) € 

^(oo)(^). 

The points t^, j = 1, . . . , fcj, represent singularity points of the functions 
at hand (note that / may be discontinuous at tj) and our main purpose is 
to establish a calculus adapted to such functions. 

Our basic hypothesis is the following. 

Hypothesis 2.2. For every z e N, the conditional law of Vi with re- 
spect to Qi is absolutely continuous on {ai,bi) with respect to the Lebesgue 
measure. This means that there exists a Qi x ;S(M)-measurable function 
Pi =pi{uj,x) such that 



E(ev(F^)i(a„fe,)(F^)) = E(^e y^^(x)p,(w,x)i(,^,fe^)(x)(ix 

for every positive, ^j-measurable random variable and every positive, 
measurable function ^ : M ^ M. 

We assume that pi G Ci{Bi) and dy\npi{uj,y) € L(^^-^{A). 

In concrete problems, we consider random variables Vi with conditional 
densities pi and we then take t], i = 0, . . . , /cj+i, to be the singularities of pj. 
This means that we choose Bi in such a way that pi satisfies Hypothesis 2.2 
on Bi. This is the significance of Bi (in the case where pi is smooth on the 
whole M, we may choose Bi = R). 

For each i € N, we consider a Qi x ;B(M)-measurable and positive function 
TTj : X M — > such that Tri{uj,y) = for y ^ {ai,bi) and TTj S Ci{Bi). We 
assume the following. 

Hypothesis 2.3. ttj g L(^^){A) and vr- G L(i+)(74). 
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These will be the weights used in our calculus. In the standard Malliavin 
calculus, they appear as renormalization constants. On the other hand, pi 
may have discontinuities at t^, j = 1, . . . , fcj, and this will produce some bor- 
der terms in the integration by parts formula; see (2.2). We may choose vTj 
in order to cancel these border terms (as well as the border terms in Oj and 
bi). 

2.2. Differential operators. In this section, we introduce the differential 
operators which represent the analogs of the Malliavin derivative and of the 
Skorohod integral. 

Simple functionals. A random variable F is called a simple functional if 
there exists some n and some Q x ;B(M")-measurable function / : 17 x M" R 
such that 

F = f{cv,Vi,...,Vn). 

We denote by S(^n,k) the space of simple functionals such that / G Cn,k, and 
5(„ fc)(j4) will denote the space of simple functionals such that / G Cn,k{^)- 
We will use the notation dy^F := ■§§-{^, Vi, ■ ■ ■ , Vn), i = 1, . . . ,n. 

Simple processes. A simple process of length n is a sequence of random 
variables U = {Ui)i<n such that 

Ui{uj) = Ui{u}, Vi(u;), . . . , y„(u;)), 

where Uj : 17 x M" ^ M, i S N, are Q x ;B(M")-measurable functions. We denote 
by 'P(n,k) [resp., P(„ fc)(j4)] the space of simple processes of length n such 
that Ui € Cn^k, i = I, . . . ,n [resp., Uj G Cn,k{A), i = 1, . . . ,n]. Note that if 
U G V(n,k), then Ui G S(^n,k) and if f/ G V(^n,k){A), then Ui G 5(„_fc)(A). 
On the space of simple processes, we consider the scalar product 

n 

{U,V)^ ■.= Y,7T,{u;,V,)Uiiij)V,iij). 
1=1 

We now define the differential operators which appear in Malliavin's calcu- 
lus. 

The Malliavin derivative. D : — > 'P(n.o)- if = /(^^i • • • i then 
DiF:=^{u;,ViiLo),...,Vniu;))lB,^^){Vi), 
DF ={DiF),<^er^n,o)- 
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The Malliavin covariance matrix. Given F = (F^ , . . . , F'^), F^ = 
f^{u>, Vi, . . . , Vn) G 'S(n,i)i the Malliavin covariance matrix is 

n 

a% = {DF\ DF^)^ = ^ 7rp(a;, Vp) dpf dpfico, Vi, . . . , K). 
p=i 

This is a symmetric, positive-definite matrix. 

The Skorohod integral. We define 5:^(^,1) <S[nfl)'- for U = (f/i)i<i<n 
such that Ui{uj) = Ui{oj, Vi, . . . , V^), we define 

6i{U) •= ~(^^(^*^*) + {T^iUi) d\np^{uj,Vi, . . . ,Vn), 

i=l 

The border term operator. For F = f{uJ,Vi, . . . ,Vn) G <5(n,o) U = 

{ui{uj,Vi, . . . ,Vn))i=i,...,n ^'P{n,o), we define 

n 

[F, U]^ = J2^Mf X Ui) X TTj X Pi) 

j=l 

n ki 

= T.Jl((f^ ^1' • • • ' Vj-l,t{-,Vj+i, Vn){7T,Pi){u;, tl-) 

1=1 j=l 

-{fx Ui)iuj,Vi,...,Vj-u 

t{ + ,Vj+l,...,Vn){^T^Pi)iiO,ti + )) 

n 

+ ^if X Ui){uJ,Vi, . . .,Vj-i,bi-,Vj+i,.. .,Vn)ilTiPi){uj,bi-) 

1=1 

n 

- ^{f X Ui){uj,Vi, . . . ,Vj^i,ai+,Vj+i, . ..,Vn)iiTiPi)iuj,ai+). 
1=1 

Remark 2.1. If we choose ttj such that 

7ri{uj,tl+) = TTi{u;,tl-) =0, i=l,...,n,j = l,...,ki, 

(2.3) 

TTi{uj,ai+) = TTi{uj,bi-) = 0, i = l,...,n, 

then [F, [7],^ = for every F G S(^n.i) ^iid U E 'P{n,i)- So there will be no bor- 
der terms in the duality formula or in the integration by parts formula. This 
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is one possible reason for introducing the weights. The other one concerns 
renormahzation. 

In our framework the duahty between 6 and D is given by the fohowing 
proposition. 

Proposition 2.1. Let F e S(^n,i) and U e Suppose that for every 

i = l,...,n, 



(2.4) E{\F 6iiU)\l a) +E{Tr,{u;, V,)\DiF x U,\1a) < oo. 
Then E{\[F,U]^\1 a) <oo and 

(2.5) E{{DF, U)Aa) = E{F6{U)1a) + E([F, U]Aa). 
If (2.3) holds true, then 

E{{DF,U)Aa) = E{F6{U)1a). 

Proof. Since vTj = on {ai^biY, we have 
E{{DF,U)^1a) 

= e(^E{^,{uj,V,)D,F X U, I G,)l}^ 

= e(i^^ / {T^iUidif){ui,Vi, . . . ,Vi-.i,y,Vi+i, . . . ,Vn)Pi{uj,y) dy 



Using integration by parts [see (2.2)], we obtain 

dif X (vTiUi) X Pi 



k. 



hi „ 

Ti{Ii{f X Ui)'KiPi) - V / . . / X {di{-KiUi) y.pi + (TTiUi) X dpi 

Ti{Ii{f X Ui)7riPi) - fx {di{-KiUi) + -KiUidlnpi) x pi. 



By (2.4), we have 

/ {\uidif\TTiPi){u:,Vi, . . .,Vi-i,y,Vi+i, . . . ,Vn) dy < oo, 

JR 

{\f{di{7riUi) + TTjUj (91npi)| x pi){uj, 14, ... , Vi-i,y, Vi+i, . . . ,Vn) dy < oo, 
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for almost all to (z A. Thus, the above integrals make sense. Since Ti[Ii[f x 
Ui)T^iPi) is the sum of these two integrals, we also obtain E(|rj(/j(/ x Ui)TTiPi)\ 
1a) < oo so that E(|[F, f/]^|lA) < oo. 

Using the definition of pi , we return to expectations and obtain 

I {-KiUi dif){uj, Vi,.. .,Vi-i,y, Vi+i, . . . , Vn)pi{uj, y) dy 

Jai 

= E{F6i{u) I g^) + r,(/,(/ X u^)7Tm). 

One sums over i and the proof is complete. □ 

Corollary 2.1. Let Q e S(^n,i){^) satisfy 

(2.6) E{lA{\7TiQ\ + \dvMQ)\)^+")<(x^, i = l,...,n, 

for some i] > 0. Then for every F G 5(„^i)(j4), U G ■P(„ i)(j4), one has 

(2.7) E{Q{DF, U)Aa) = E{F5{QU)1a) + E([F, QU]^1a). 

Proof. We need only check that F and U = QU satisfy (2.4). We have 

\6i{QU)\ < \dvMQ)\m + \7riQ\{\dvm + \Ui\\dlnp,\). 

Since U G P(„^i)(A), one has Ui^dvJJi G L(^^-^{A) and by Hypothesis 2.2, 
d\npi G L(oo)(^). So, using (2.6), we have 6i{QU) G L(^i^-^{A) and since F G 
L(^)(A), we obtain E{F6i{QU)\) < oo. 

We have DiF, Ui G L(oo)(^) and vrjQ G L(i+)(A), so E(7rj|A-F' x (QC/i)|) < 
oo. □ 

The Ornstein-Uhlenbeck operator. We now introduce L := (5(-D) : '5(n,2) ^ 

'5(n,0)- 

n 

LF := - ^(9.(7r, 9,/) + vr, difdlnp,)iu;, ^i, . . . , K) 

i=l 

n 

= -E(^^+^^^ln^'05J + ^*52/)(a;,yi,...,K). 

As an immediate consequence of the duality relation (2.5), we obtain the 
following. 

Lemma 2.1. Let F,G £ '5(n,2) A £ Q. Suppose that for every i = 
l,...,n, 

E[{\FLiG\ + \GL,F\ + mlDiF x D,G\)1a] < oo. 
T/ien E(|[F,Z)G]^|1a) < oo, E(|[G, L>F]^|1a) < oo and 
E{FLG1a) + E{[F,DG]Aa) =E{{DF,DG)Aa) 

= E{GLF1a) + E{[G,DF]^1a). 
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We denote by C^(M'^) the space of the functions : M*^ — > M which are k- 
times differentiable and such that (j) and its derivatives of order less then or 
equal to k have polynomial growth. Standard differential calculus gives the 
following chain rules. 

Lemma 2.2. (i) Letcpeel{W^) and F = {F\...,F'^), E 5(„_i)(^). 
Then (l){F) G S^ri,i){A) and 

d 

(2.8) D(l){F) = Y,dk^{F)DF''. 

k=l 

(ii) If cl) £ C^iW^) and F' eSi^n,2){A), then 4>{F) e S(n,2){A) and 

k=l k,p=l 

(iii) Let F e cS(„_i)(A) and U £ Then FU G Vi^n,i)i^) and 

5{FU) = F5{U) - {DF,U)^. 

In particular, if F € S(^n,i)i^) and G £ S(^n,2){^) > then FDG GV(^n,i)i^) and 

(2.9) 5{FDG) = FLG- {DF,DG)^. 

Remark 2.2. Let us define Ll^^{A) to be the closure of "^(^,0) "^ith 
respect to the norm associated with the scalar product {U, V) = E((?7, V).,^). 
If [F, U]t, is not null, then the operator D : 5(ri,i) C L^(r2) ^(n,o) -^w,n(^) 
is not closable. 

Suppose, for example, that Vi is exponentially distributed and Vi, i = 
2,...,n, are arbitrary and independent of Vi. We take vri = 1 and tTj = 
0, i = 2, . . . ,n. So we define our calculus with respect to Vi only. In this 
case, ai = 0, 6i = oo and there are no points tj . Now take Fn = fniVi) with 
fn{x) = 1 — nx for < x < 1/n and fni^) = for x > 1/n. Also, take ui{x) = 
1 — X for < X < 1 and ui{x) = for x >1 and write the duality formula 
E((DF„, U)^) = E{Fn5iU)) + E{[Fn, U]^). Since [F„„ f/], = 1 and F„, ^ in 
L'^i^l), we obtain lim„^oo E((i:>F„, C/)^) = 1 and so DFn in Ll^^{A). 
This proves that D is not closable. 

But if C/Jtt = for every F, U [this happens, e.g., if we choose vTj to 
satisfy (2.3)], then the duality formula (2.5) guarantees that D and 6 are 
closable. But we will remain at the level of simple functionals and will not 
discuss the extension to the infinite-dimensional setting. 

Remark 2.3. The above differential operators and the duality formula 
(2.5) represent abstract versions of the operators introduced in Malliavin cal- 
culus and of the duality formula used there. In order to see this, we consider 
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the simple example of the Euler scheme for a diffusion process, correspond- 
ing to the time grid = sq < si < ■ ■ ■ < Sn = s. This is a simple functional 
depending on the increments of the Brownian motion B, that is, Vi = B{si) — 
B(si-i), i = 1, . . . ,n. The variables on which the calculus is based are inde- 
pendent Gaussian variables. It follows that pi{uj,y) = (27r(sj — Si_i))^^/^ x 
exp(— ?/^/2(si — Si-i)). Since pi is smooth on the whole of M and has null limit 
at infinity, there will be no border terms, so we take Oj = —oo, bi = oo and 
= O.U F = f{Vi, . . . , Vn), then DiF = dif{Vi, . . . , K) = :D,F1[,^_^,,^)(s), 
where DgF is the standard Malliavin derivative. We take vr, = — Si-i so 
that 

{DF, DG)^ = f] TTiD.FDiG = D^FD^G du. 

We note that here the weights are used in order to obtain the Lebesgue 
measure. Moreover, we have dy\npi{y) = —y/{si — Sj_i), so 

n 

= -Y.{diUi{Vi,. . .,Vn){s^ - Si-i) - Ui{Vu . . . , 
1=1 

So we find out the standard Malliavin calculus. 

Remark 2.4. If [F, GJ^r = 0, the calculus presented here fits into the 
framework introduced by Bouleau in [5]: in the notation there, the bilinear 
form (F, G) {DF, DG)t^ leads to an error structure. A variety of examples 
and applications of these structures is discussed. That framework mainly 
focuses on the error calculus, but examples of applications to sensitivity 
computations are also given and an integration by parts formula is derived. 
This works well in the particular case of a one-dimensional functional. More- 
over, the differential calculus is based on a single noise Vi as in Corollary 2.2 
below (so the weights vTj do not come into the nondegeneracy condition). In 
a more general framework, the nondegeracy condition involves the weights 
TTj, i G N, and a more detailed analysis must be undertaken (see the following 
section). 

2.3. The integration by parts formula. We consider F = {F^, . . . , F'^) G 
^{n 1) (^) ^^'^ define 

SpiA) ■.= {G = aF xQ:QeSf^^^)iA),Qi satisfy (2.6)}. 

We think of G G @f{A) as a random direction in which F is nondegenerate 
(in Malliavin's sense). 

The basic integration by parts formula is given in the following theorem. 
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Theorem 2.1. Let F = {F\ . . . , F"^) £ Sf^^^^iA) and G eOpiA), G = 
(Tp ^ Q- Then 



sij^Q'DF' 



\i=l 



^{F),Y,QWF^ 



i=l 



€i(l+)(^) 



and for every cj) G C^iW^), one has 

E{{Vct^{F),G)lA) = E((l>{F)5(Y,Q'DF'\l, 



(2.10) 



i=l 



+ E 



(t>{F),Y,QWF' 



1=1 



Proof. Using (2.^ 



{D<p{F),DF\ = Y,d,cl>{F){DF^,Dr)^ = J2dj^{F)a'i 



Since G = ap x Q, we obtain 



{y<p{F),G) = y: dj<piF)G^ = y: d,m y: q'4 

i=i i=i 



i=l 



= jZQ'Y dMF)a'^ = Y Q\D(p{F),DF%. 

i=l j=l 1=1 

Note that <j){F) G and € 7^(„_i)(A). Since the Qi satisfy (2.6), 

one may use the duahty formula (2.7) and thereby obtain (2.10). □ 

We now give a nondegeneracy condition on cr_p which guarantees that all 
of the directions are nondegenerate for F. 

We assume that det ap on A and define 7_f = ap^. We also assume that 
7r/(det 7_p)^, VT; det 7_f, 7r/7r;(det 7i;')^ G L(i_^)(j4) for every I = 1, . . . ,n. This means 
that there exists r/ > such that 



(2.11) E[lA(|vr;|(det7F)^ + k;i(det7i. + |7r;|(det7F)2))^+''] < oo. 



Lemma 2.3. Assume that (2.11) holds true and that F G S^^ 2)(^)- Then 
Sf^,^{A)CQp{A). 
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Proof. Let G G Then G = ap x Q with Q = -fp xG. We 

write = ap x det7i;', where ap is the algebraic complement. It follows 
that Q' = det7F x S' with S' = J^Ui ^^^F- 

Let us check that (2.6) holds true for Q^, i = 1, . . . ,d. Since tt; S L(^^^{A) 
and DiF^ £ L(^f^-^{A), one has ap, deters G L(^^^(A) and since G-^ € L(^^-^{A), 
we have 5* G L(oo)(74). Moreover, by (2.11), 7r/det7F G L(i+)(A), so vr^Q* = 
(7ridet7i.)5'GL(i+)(A). 

We now check that Di{-kiQ^) G L(i_|_)(^). We write 

n 

Dia% = tt'iDiF'DiF^ + ^ T^kDiiD^F' D^F^). 

k=l 

Since F G S('^„^2)(^)' ^e have DiFWiF\ DiiD^FW^F^) G L(oo)(^) and 
consequently Dia^p = 9i + 621t'i with 6i, 62 £ L^^-j{A). Then Z)i(deto"i7) = 
/i + i^7r[ and DiS'' = /ij + fivr^ with /x, u, ^i, Vi G L(oo)(^). 

Using (2.11), we obtain 

Di{-KiQ')=-K[dei-iFS'-TTi{det-iFfDi{deiaF)S' + Tridei-iFDiS' 

= 7r;det7j7'5* — 7r;(det7^)^(;U + i^vr^)5* + vr^ det7i?(/Xj + fjvr^) 

ei{i+)(A) 

and the proof is complete. □ 

As a consequence, we obtain the following. 

Theorem 2.2. Let F = {F^ , . . . ,F'^) £ Sf^^^^^A) and G G 5(„^i)(^). Sup- 
pose that (2.11) holds true. Then 



6iGY.j^;DF^ 
\ j=i 



ct>{F),GY,j^;DF^ 



e^{i+)(^) 



and for every <p G Cp(IR'^), one has 



E{di(l){F)GlA) = E 



(t>{F)6(Gj2l';DFAl, 



+ E cP{F),GY,-fj;DF^ 

V L i=i 

/or e?;erj/ i = 1, . . . , d. 

Suppose that iri, I = 1, . . . ,n, satisfy (2.3). Then 

E{dict>{F)GlA) = E{cl){F)Hi{F, G)1a) 
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with 



H,{F, G) = 6[ Gj2-f^;DF^ = Y^iG^'^LF^ - {D{Gj^;),DF^)^). 



j=i J j=i 



Proof. We take G= (0, . . . ,0, G,0, . . . ,0) with G occupying the ith 
place so that di(p{F)G = (V0(F),(5). In view of Lemma 2.3, G G QpiA) 
and G = ap x Q, with = G'fp. One then employs Theorem 2.1. In order 
to obtain the second equality in the expression for Hi {F, G) , one employs 
(2.9). □ 

There is one particular situation in which the nondegeneracy condition 
(2.11) does not involve the weights — when if F is one-dimensional and the 
integration by parts formula is based on a single random variable Vi. We 
then have the following corollary. 

Corollary 2.2. Let F = f{Vi,. . . ,Vn) e S(n,2){A) anc? G G cS(„,i) (A). 
Suppose that there exists some I G {1, . . . ,n} such that 



for some 7] > 0. Consider the weights vTi = for i ^ I and let tt/ be an 
arbitrary function which satisfies tt/ G Lj^^^{A) and vr^' G L(i_|_)(^). Then 



5{GjfDF), [(j){F),G'yFDF]^ e L^i_^){A) and for every cp € Cl{R) , one has 
(2.13) E(0'(F)G1a) = E{(l)iF)5{G^FDF)lA) + E([0(F), G-fFDF]AA). 

Proof. Note that cjf = 'n-i{Vi)\DiF\'^ . We return to the proof of Theo- 
rem 2.1 and write G = Qgf with 



Then i^iiYijQ = ^(Vi, . . . , Vn) /\DiF\'^ and, as a consequence of the hypothesis 
(2.12), one has 'Kiiyi)Q, 9y.(7r(V;)(5) G L(ij^-^{A), i = 1, . . . ,n. So we may use 
the duality relation to conclude the proof. □ 

3. Pure jump diffusions. In this section, we will use the integration by 
parts formula presented in the previous section for a pure jump diffusion 
iSt)t>o- We will use the notation from [14]. We consider a Poisson point 
measure N{dt, da) on M with positive and finite intensity measure ^{da) x dt, 
that is, E(A^([0, t] x A)) = n{A)t. We denote by Jt the counting process, 
that is, Jt ■= N{[0,t] x M) and we denote by Tj, i G N, the jump times of Jf. 



(2.12) 



E[1a(AF)~'^'+''^] < oo 




if 7riiVi)\DiF\^^0, 
if TTiiVi)\DiF\^ = 0. 
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We represent the above Poisson point measure by means of a sequence Aj, 
« E N, of independent random variables with law u{da) = x ^[da). 

This means that iV([0, t]x A) = card{Ti < t : € A}. 
We look at the solution St of the equation 



We will work under the following hypothesis. 

Hypothesis 3.1. The functions ^ c(t,a,x) and x ^ g{t,x) are 

twice differentiable and have bounded derivatives of first and second or- 
der. Moreover, we assume that they have linear growth with respect to x, 
uniformly with respect to t and a, that is, \c{t, a,x)\ + \g{t, x) \ < K{1 + 

On each set {Jt = n}, St is a simple functional of Ai, . . . , A„ and Ti, . . . , T„. 
In the first subsection, we present the deterministic calculus which permits 
us to compute the Malliavin derivatives and in the following two subsections, 
we give the integration by parts formula with respect to the amplitude of the 
jumps and with respect to the jump times, separately. Finally, in the fourth 
subsection, we briefly present the mixed calculus with respect to both. 

We shall remain in the one-dimensional case (i.e., St S M) because the 
multi-dimensional case is more involved from a technical point of view. Our 
purpose is to illustrate the way in which the integration by parts formula 
works for Poisson point measures and to emphasize the specific difficulties. 
The heavy techniques related to the multi-dimensional case would obscure 
these specific points, but the machinery works just as well in this case. 

3.1. The deterministic equation. We fix some deterministic = uq < 
ui < ■ • • <Un <T and define u = (ui , . . . , n„) . We also fix a = (ai , . . . , a„) G 
M". To these fixed numbers, we associate the deterministic equation 



where Jt{u) = k if Uk < t < Ufc+i. We denote by st{u,a), or simply by st, 
the solution of this equation. This is the deterministic counterpart of our 
stochastic equation. On the set {Jt = n}, the solution St of (3.1) is repre- 
sented as St = st{Ti, ... ,Tn, Ai, An) . 




(3.2) 
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In order to solve this equation, we introduce the flow $ = ^uii,x), < 
u<t, X S M, which solves the ordinary integral equation 

(3.3) ^u{t,x)=x+ f g{r,^u{r,x))dr, t>u. 

J u 

The solution s of the equation (3.2) is given by 
So = X, 

(3.4) St = {t, ) for Ui<t< Ui+i , 
Su,+, =s^' +c{ui+i,ai+i,s-) 

i-f-1 i+l 

= ^Ui{ui+i,SuJ +c(nj+i,aj+i,$„.('Uj+i,s„.)). 

Our aim is to compute the derivatives of s with respect to uj, aj, j = 1, . . . ,n. 
We first introduce some notation. We define 

eu,t{^) :=exp^y dxg{r,^u{r,x)) drj . 

Since ^mif, s^.) = Sr for Ui <r < Uj+i, we have 

eui,t{su,) = exp (^J dxg{r, Sr) drj for m < t 



< Ui+l- 



Since 

da:^u(t,x) = 1+ I da:g{r,^u{r,x))dx^u{r,x)dr, 

Ju 

it follows that 

dx^u{t,x) = eu,t{x) 

and since 

du^u(.t,x) = -g{u,x) + I dxg{r,^u{r,x))du^u{r,x)dr, 

J u 

we have 

du^uit.x) = -g{u,x)eu,t{x). 

Finally, we define 

q{t,a,x) := {dtc + gdxc){t,a,x) +g{t,x) - g{t,x + c{t,a,x)). 
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Lemma 3.1. Suppose that Hypothesis 3.1 holds true. Then st{u,a) is 
twice differentiahle with respect to Uj,j = 1, . . . ,n, and with respect to aj,j = 
1,. . . ,n, and we have the following explicit expressions for the derivatives. 

A. Derivatives with respect to Uj. For t <Uj, dujSt{u,a) = 0. Moreover, 

du.Suj = idtc + g{l + dxc)){uj,aj,Suj-). 

For Uj <t < Uj+i, 

dujSt = q{uj,aj,Suj-)euj,t{suj), 
(3.5) dujSuj^i— — Q{uj , Oj , Suj—)euj ^Uj^i{suj) 1 

duj ^Uj-^^i — Qi^j 1 ; — ) (1 ~l~ dxC(^Uj-\-l , + 1 ) "SMj^i — ,Mj + i i^Uj ) • 

Finally, for p> j + 1 and Up<t < Up+i , we have the recurrence relations 
(3.6) 

dujSup+i = (1 ~^ ^xc{up^l,ap-^-l, Sup^i-))eup,up+i{sup) dujSup. 
Define T{f) '■= dtf + g dxf . The second order derivatives are given by 
dljSuj^ =T{g){uj,aj,Su^-), 
^Ij^uj =T{dtc + g{l + dxc)){uj,aj,Su^-). 

Define 

Pj{t) = du^Cu^^tiSuj) 

= euj,t{su,)(y-dxg{uj,Suj) + q{uj,aj,Suj-.) J dlg{r, Sr)eu^^r{suj) dr^ . 

Then for Uj <t < uj^i, 

dl,st{u, a) = T{q){uj,aj,Su^-{u, a))eu^^t{suj ) + q{uj,aj,Suj-{u, a))pj{t) 
and 

dl^ Su,+^ = T{q) {uj , Oj , s„^,_ ) (1 + dxc) {uj+i , aj+i , s„ .+,_)e„ . (s„ . ) 
+ q'^{uj,aj,Su^-) dlc{uj+i 

+ Qiuj , aj ,Suj-~)il + dxc) {uj+i , flj+i , )pj (uj). 

For p>j + l, we define 

Pj,p{t) = du^eup,t{sup) = eup,t{sup)du,Sup f dlg{r, Sr)eup,r{sup) dr. 

J Up 



INTEGRATION BY PARTS FORMULA FOR LOCALLY SMOOTH LAWS 17 

Then for p> j and Up <t < Up^i, we have the recurrence relations 
dljSt = eup,t{sup) dl.Sup + pj,p{t,u,a) du^Sup, 

^UjSup+i = d^c{Upj^i, ttp+i, Sup^i-){eup,Up+i{Sup) dujSup) 

X {pj,p{Up-\-l) dujSup + Gup,Up+i{Sup) d^.Sup) ■ 

B. Derivatives with respect to aj . For t <Uj, da^ {u, a) = and for t>Uj, 
dajSt{u,a) satisfies the equation 

Jt{u) 

daj St = daC{ Uj , CLj , Suj — ) ~l~ •^Ui — 

i=j+l 

(3.7) 

+ / dxg{r,Sr)da^Srdr. 

Juj 

The second-order derivatives satisfy 

Mu) 

dljSt = dlc{uj,aj,Su^-) + E dlc{ui,ai,Su,-){dajSui-f 

i=j+i 



(3.8) + / dlg{r,sr){da^srfdr 

Mu) . 

UiiChii Sin — 

) dlj Su,- + / da:g{r, Sr) dl^ Sr dr. 

i=j+l 

Proof. A. It is clear that for t < Uj, St does not depend on Uj and so 
dujSt = 0. We now compute 

dujSu^- =du^^Uj^j{Uj,Su^^j) =g{uj,^uj.Auj,Suj.i)) =g{Uj,Suj-). 

Then 

duj Suj — {Suj — ~l~ c{^Uj , CLj , — ) ) 

= dtc{uj,aj,Suj-) + {l + dxc{uj,aj,Su^-))du^Su^- 

= dtc{uj,aj,Suj-) + {1 + dxc{uj,aj,Suj-))g{uj,Suj-). 
For Uj < t < Uj^i, we have 

dujSt = du^^uj(t,suj) = euj^t{suj){-g{uj,Suj) + 9„^s«,) 
= euj,t{suj){-g{uj,Suj) + dtc{uj,aj,Su^-) 

+ {l + dxc{uj,aj,Suj-))g{uj,Suj-)) 
= euj,tisuj)qiuj,aj,Suj-) 



18 V. BALLY, M.-P. BAVOUZET AND M. MESSAOUD 

and the same computation gives dujSuj+j^- = euj,Uj+i{suj)qiuj,aj,Su^-). Fi- 
nally, 

dujSuj+^ = {l + dxc{u 

= {l + dxc{uj+i,aj+i,Suj^j^-))eu^^Uj+iisuj)q{uj,aj,Suj-). 
We now assume that Up <t < Up+i, p > j + 1, and we write 

duj St — ) ) — ^Up,t i^Up ) Sup • 

The same computation gives du^Sup+i- = eup,up+i{sup) du^Sup. Finally, we 
have 

dujSup — diij{Sup— c{^Up, dp, Slip—)) — (1 ~\~ dxCi^Up, CLp, Sup^}^ S,i^^ — 
= (1 + dxC{Up, Op, Sup — y)^Up-^,Up (•S-Up-i ) duj Sup-i 

and the proof is complete for the first-order derivatives. The relations con- 
cerning the second-order derivatives are obtained by direct computations. 

B. Using the recurrence relations (3.4), one verifies that for every t G [0,T], 
Uj st{u,a) is continuously differentiable and one may then differentiate in 
equation (3.2) (this was not possible in the case of the derivatives with 
respect to uj because these derivatives are not continuous). □ 

As an immediate consequence of the above lemma we obtain: 

Corollary 3.1. Suppose that Hypothesis 3.1 holds true and that the 
starting point x satisfies \x\ < K for some K . Then for each n € N and 
T > 0, there exists a constant Cn{K,T) such that for every < ui < ■ ■ ■ < 
Un<T, a € M" andO<t<T, 

(3.9) max {\st\ + \du^st\ + \dlst\ + \da,st\ + \dlst\)iu,a) < Cn{K,T). 

j=l,...,n J J J J 

Finally, we present a corollary which is useful in order to control the 
nondegeneracy. 

Corollary 3.2. Assume that Hypothesis 3.1 holds true and there exists 
a constant i] > such that for every {t,a,x) G [0,T] x M x M, one has 

\l + dxc{t,a,x) \ > rj, 

(3.10) 

\q{t,a,x)\ > T]. 

Let n £N be fixed. Then there exists a constant e„, > such that for every 
j = 1, . . . ,n and every (n, a) G [0, T]" x M", we have, 

(3.11) inf \du,st{u,a) \ >e„. 

t>Uj ^ 

Proof. Since dxQ is bounded, there exists a constant C > such that 
es,t{x) > e"*^^ for < s < t < T. One then employs (3.5) and (3.6). □ 
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3.2. Integration by parts with respect to the amplitudes of the jumps. In 
this section, we wih use the integration by parts formula for St which will 
be regarded as a simple functional of Aj,z € N. So, with the notation from 
Section 2, we have Vi = Aj and Q = a{Ti :i S N}. We recall that Jt = n on 
{Tn < i < Tn+i}- Then, on { Jt = n}, we have 

St = st{Ti, ...,Tn,Ai,.. . , A„), 

where st is defined in the previous section [see (3.2)]. 

We assume that Hypotheses 3.1 and 2.1 [i.e., E(|Aj|P) < cxd for all p G N] 
hold true. Moreover, we consider some qo < qi < ■ ■ ■ < qk+i and define 

k 
1=0 

We assume the following. 

Hypothesis 3.2. The law of Aj is absolutely continuous on / with 
respect to the Lebesgue measure and has the density p{y) = lj{y)e/^y\ that 
is, 

E(/(A,)) = |^/(y)e''(J^) dy 

for every measurable and positive function /. 

The function p is assumed to be continuously differentiable and bounded 
on /. 



Therefore, Hypothesis 2.2 holds true. 

Since p is not differentiable on the whole of M, we work with the following 
weight. We take a € (0, 1) and f3 > a and we define 

{qi+i -yTiy -QiT, for ye{qi,qi+i),i = 0,...,k, 
0, for y G (goWfc+i)''- 



Tr{y) 



We introduce the following convention: if 5 = qk+i = +oo or a = qo = — oo, 
we define 



7r(y) 



{y-qk)"\y\ ^, for y 

{qi -yT\y\~'^, fory<gi. 

Since a G (0,1), we can show by elementary computations that vr satisfies 
Hypothesis 2.3, that is, vr € L(^^^{A) and vr' € L(i4.)(A). 
Let A := {Jt = n}. In view of (3.9), 

(ai, . ..,«„,) ^ stiTi{uj), . . . ,r„(u;),ai, . . . ,a„) 

is twice continuously differentiable and has bounded derivatives. It follows 
that St G S'(„^2)(^)- 
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The differential operators which appear in the integration by parts for- 
mula are 

DiSt = daiStiTi, . . .,r„, Ai, . . . , A„), 

n . 

LSt = -Y.U{A,)dlst{Ti,...,Tn,Ai,...,An) 

1=1 ^ 

+ (^tt' + 7r^yAi)daMTi,---,Tn,Ai,...,An)^, 

n n 

as, = ^7r(A0| A5f P = ^7r(A,)|9a,st(ri, . . . , r„„ Ai, . . . , AJp, 

i=l 4=1 
1 1 



as, Er=i 7r{Ai)\d,MTi, . . . , T„, Ai, . . . , A„)|2 " 

All of these quantities may be computed using (3.7) and (3.8). 

The result which is used in sensitivity computations is the following. 

Proposition 3.1. Suppose that Hypotheses 3.1 and 3.2 hold true and, 
moreover, suppose that there exists a positive constant rj such that for every 
t,a,x, we have 

(i) \dac{t,a,x)\>r], 

(3.12) 

(n) \l + dxc{t,a,x)\>r]. 

Take a S (0, 1 /2) and (3> a. Then for every dijjerentiable function : M — > R 
which has linear growth and for every n>l, 

E{4>'{St) 5,Stl|j,=„|) = E{4>{St)Hnl{j,=n}), 
with, on {Jt = n}, 

Hn '■= Hn{St,dxSt) 

= dxSas,LSt-^sADSt,D{dxSt))^ - dxSt{DSt,D^s,).- 

Proof. Let n G N* be fixed. We already know that St G 5(„ 2)(^) with 
A = {Jt = n}. 

Moreover, dxSt = dxSt(Ti, . . . , T„, Ai, . . . , A„) and dxSt is computed by 
the recurrence relations 

dxSo = 1, 

dxSt = {l + dxc{ui,ai,Su,--))dxSu,~+ f dxg{r,Sr)dxSrdr, Ui<t<Ui+i. 

Jui 
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It is then easy to check that dxSt and its derivatives with respect to aj, 
i = 1, . . . , n, are bounded on the set {Jt = n} and consequently that dxSt € 

S(^n,lM). 

• Suppose that n = 1. We will use Corollary 2.2, so we check that the 
nondegeneracy condition (2.12) holds true. One has 



da^St = daC{ui,ai,Sui-) + I dxg{r,Sr)da^Sr 

J Ml 

so that, using (3.12), 

\daiSt\ = \dac{ui,ai,Su^-)\ex:p( f dxg{r,Sr) 



dr 



> c 



for some positive constant c. Inequality (2.12) then follows. Then the in- 
tegration by parts formula (2.13) holds true for St and dxSt on the event 
^ = {Jt = !}• Moreover, by our choice of vr, the border terms are canceled, 
which gives E((/)'(5t)a^5a{j,=i}) =E(0(5t)lfil{j,=i}) with 

On {Jf = 1}, we have 

7r{A^)dxSt{Ti,Ai)daMTu^i) 



7r{Ai)dxSas,DSt 



7r(Ai)|a,,Si(ri,Ai)|2 
dxSt{Ti,Ai) ^ 

St{Ti,Ai) 



• Now suppose that n > 2. In this case, we will use Theorem 2.2, so we look 
at the nondegeneracy condition (2.11). Since vr is bounded, this amounts to 
finding 6 > such that for i = 1, . . . , n, 

(3.13) m{J.=n}{{l + |vr'(A,)|)7lj'+'] < oo. 

We recall that Tr{y) = ^^Lgfe+i " vTiv - ^i)"lfe,g,+i)(2/)> so 

10, ifyG(go,9fc+i)^ 

We choose 5 > Q such that 2a(l + 5) <l and (1 - a)(l + 5) <l [which is 
possible because a G (0, 1/2)]. In particular, since p is bounded on / and Aj 
have finite moments of any order, this gives 

E(7r(Ai)-2(i+<5))<oo and E(|7r'(Ai)|i+'^) < oo. 

The proof of (3.13) is different for i = n and i = 1, . . . , n — 1. 
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First, take i < n. One has \da„st \ = |5ac(n„, exp(/^^ dxg{r,Sr)) > c, 

so ast > c^7r^(A n)- Since Aj and A„ are independent, 

E[1|^,=„}((1 + |7r'(A,)|)7lj'+'] < c-'E[l{j^^^y{{l + |7r'(A,)|)7r-2(A„))^+^] 

= c-2E(7r-2(i+^)(A„))E[(l + |7r'(Ai)|)'+'] 
< oo. 

Now, take i = n and write ast > ''r^(A„_i)|L'„_i5tp. A simple computation 
shows that 

da„_iSt = dac{un-i,an-i,s - ) (1 + 9^c(n„, a„, s - )) cxp f / dxg{r,Sr)dr 

Using (3.12), we obtain da„_iSt > c> 0, so ast > c^7r^(A„_i). Consequently, 

E[l|,,=„}((l + |^'(A0|)7lj^+'] 

< c-2e[1|^,=„|((1 + |7r'(A,)|)7r-2(A„))^+^] 

= c-2E(7r-2(i+'5)(A„_i))E[(l + |7r'(A„)|)i+^] < oo 
and the proof is complete. □ 

Remark 3.1. Suppose now that p is difFerentiable on the whole of M. 
We then take no weight, vr = 1 and hypothesis (3.12)(i) gives as^ > c on 
{Jt = n} for all n € N*. So the above integrability problems disappear. In 
particular, hypothesis (3.12)(ii) is no longer necessary. This case is discussed 
in [1]. 

3.3. Integration by parts with respect to the jump times. In this section, 
we differentiate with respect to the jump times Tj, i G N. It is well known 
(see [2]) that conditionally to {Jt = n}, the law of the vector (Ti, . . . ,T„) 
is absolutely continuous with respect to the Lebesgue measure and has the 
following density: 

p{uJ,ti, ...,tn) = — l{0<ii<-<t„<i}(*l;- ■ • )^n)l{Jt(a;)=n}• 
In particular, for a given i = l,...,n, conditionally to {Jt = n} and to 
{Tj,jy^i}, Ti is uniformly distributed on [rj_i(a;),Tj+i(a;)]. Therefore, it 
has the density (with the convention Tq = 0, T^+i = t) 

Pii^^u) = j,,^^^^^ ]_ l[T,-i(a;),T.+iH] (u) du, i = l,...,n. 

Since pi is not differentiable with respect to u, we must use the following 
weights: 

7ri{uj,u) = {Ti+iiuj) -u)°'{u-Ti^i{uj))°'liT,_i(uj),Ti+i{i^)]{u), i = l,...,n, 
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with a G (0,1/2). 

In order to fit with the notation from the first section, we take Vi = Ti, 
ki = 2, t\ = Ti_i and tj = Ti+i. We have G = a{Ai,i € N) V a{Jt). We fix n 
and work on the set A = {Jt = n}. Hypotheses 2.1, 2.2 and 2.3 then hold true 
and St = st{Ti, ... ,Tn, Ai{uj), An{Ld)). So St is a simple functional and 
the function which represents St is twice differentiable and has continuous 
derivatives on the whole of M". The differential operators are 

DiSt = duiSt{Ti, . . . ,r„, Ai(6j), . . . , A„(a;)), 

n 

crst = ^T^i{'^,Ti)\du^st{Ti, ...,Tn, . . . , A„(tj))p, 

1=1 

LiSt = -{ir'idu.st + vTi dl^st){Ti, . . . , T„, Ai(a;), . . . , A„(tj)) 
and all of these quantities may be computed using Lemma 3.1. 



Proposition 3.2. Suppose that Hypothesis 3.1 holds true. Suppose, 
moreover, that (3.10) is satisfied, that is, that 

\q{t,a,x)\ >r]>0, 

\{l + d^c){t,a,x)\>r]>0 

for some rj > 0. Take a G (0, ^) . Then for every n > 4 and every continuously 
differentiable function (p which has linear growth, we have 

E{(t>'{St) a,5a{j,=„}) = E{<p{St)Hnly,=n}) 
with, on {Jt = n}, 

Hn '■= Hn{St,dxSt) 

= dxSt-fs,LSt-js^{DSt,D{d.,St))^ - dxSt{DSt,D-fs,)^. 

Proof. From (3.9), we know that st{u, a) and its derivatives up to order 
two with respect to Ui, i = I, . . . ,n, are bounded on [0,T]'^. It follows that 
5iG5(„,2)(^). 

Since vTj are bounded, the nondegeneracy condition (2.11) amounts to 

E[l|,,=„|7sf ^'^] < oo and E[l{j^^^y^lf+^'>\ni{T,)\'+^>]<oo 
for some r/ > 0. 

Let us prove that E[l|j^=„|7j[^+''^ K(ri)|i+''] < oo. We define 5i = Ti- 
Tj_i and 6n+i = T — Tn, so that vTj = 5f5f_^i. We use (3.11) in order to obtain 

n n 

^5. = E S?+iS?\duMTi, . . . , Ai, . . . , A„)|2 > £2 ^ df+i^r. 

i=l 1=1 
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Since vr^(7i) = a{—6°'_^^ 6f + ), we must check that, for every i 

l,...,n, 

/ n \ -2{i+»?)- 



< oo. 



Take i = l and write 
/ n 



E 



-2(l+r,)-| 



< E[(5°-i5^)^+''(5^5^)-2(i+^')] 



: E(5("-')('+'') )E(52""^'+^^ )E(,53-'"('+'') ) . 



Since 5i is exponentially distributed of parameter ^(M), a necessary and 
sufficient condition in order to have E((5~^) < oo is p < 1. We then choose 
r] sufficiently small that (1 — q)(1 + r/) < 1 and a(l + r]) < 2a(l + ry) < 1 
(which is possible because < a < 1/2) and we have E((5[" '^)(-^+'')^ 
E(52-"(^+'') ) < oo and E(53-'"('+'') ) < oo. So 



< oo, 



E 



-2{l+7?)- 



< oo. 



We now write 
E 



-2(l+r,)-| 



= E(4"-')('+^))E(5;^('+''))E(53-'"('+^'))E(,54 



Recalling that 5i has finite moments of any order, by the choice of r], we 
obtain 

/ n \ -2{l+»?)- 

E^{S'^6r')'-'''(j2Sf^,6fj J <oo. 

Since n > 4, the same argument works for i = 2, . . . , n and leads to E[l| x 

2(l+^)l ^ 1-1 

Remark 3.2. Suppose that n = 2. Then 

/ n \ -2(l+r?) 



~a(l+r)) 



x(J- 



(a+l)(l+r?) 



2a(l+r,)^--(l-a)(l+r?)^ 
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and this quantity is not integrable for a > 0, 77 > 0. 

Remark 3.3. For n = 1, one may use Corollary 2.2 in order to obtain 
an integration by parts formula. 

But for n = 2 and n = 3, we are not able to handle the nondegeneracy 
problem. In our numerical examples, we will use the noise coming from the 
amplitudes of the jumps in order to solve the problem for n = 2 and n = 3. 

3.3.1. Examples. • We consider the geometrical model 

dSt = St{rdt + a{t,a) dN {t,a)) . 

In this case, g{t,x) = xr and c{t,a,x) = xa{t,a). It follows that 

q{t, a,x) = x dta{t, a) + xra{t, a) + xr — r{x + xa{t, a)) = x dta{t, a). 

In particular, if a does not depend on the time, the model is degenerate from 
the point of view of the jump times. The nondegeneracy condition becomes 

\dta{t,a) \ > e. 

On the other hand, the condition |1 + dxc{t,a,x)\ > rj becomes 

|l + a(t,a)| >r?. 

• We now consider the following Vasicek-type model: 

dSt = Strdt + a{t, a) dN{t, a). 

In this case, g{t,x) = xr and c{t,a,x) = a{t,a). It follows that 

q{t, a, x) = dta{t, a) + xr — r{x + a{t, a)) = dta{t, a) — ra{t, a). 

Suppose that a does not depend on the time so that dta = 0. Then the 
nondegeneracy assumption becomes 

|a(«)| > 

The condition |1 + dxc{t,a,x)\ > rj becomes 

|1 + a{a)\ > rj. 

3.4. Mixed calculus. In this section, we briefly present the differential 
calculus with respect to both noises coming from the jump amplitudes and 
from the jump times. So the random variables will be Vi = Ti,i = 1, . . . ,n, 
Vn+i = Aj,i = l,...,n, and G = (j{Jt). We combine the results from the two 
previous sections (and we keep the notation therefrom). We still assume 
Hypotheses 3.1 and 3.2. The differential operators are 
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We will use the weights defined in the previous sections, namely 

■Ki{u;,u) = {Ti+i{uj) - u)'^{u-Ti_i{uj))"l[T,_:,{uj),T,+ i(u;)]{u), i = l,...,n, 

7ri(y) =vr(y) = ^(gp+i i = n + 1, . . . ,2n, 



where a G (0, ^) 
We have 



-{'K[{Ti)du,st + -Ki{Ti)dl^St), for i = l,...,n, 

-(7r(Ai) dl^st + W + da^st), for i = n + 1, . . . , 2n. 



Finally, LSt = X)f=i LiSt- All of these quantities may be computed using the 
formulas from the previous sections. 

Theorem 3.1. Suppose that Hypotheses 3.1 and 3.2 hold true and that 

(i) |g(t,a,x)| >e>0, 

(ii) \dac{t,a,x)\>e>Q, 

(iii) |(l + a^c)(t,a,x)| >e>0. 

Then for every n>l and every continuously differentiahle function (j) which 
has linear growth, we have 

with, on {Jt = n}, 

Hn '■= Hn{St,dxSt) 

= d^Sas^LSt - 75, {DSt,D{dxSt))^ - d^StiDSt, D^s^).- 

Proof. We write 

(7 St > (vr„(a;,T„)|(9„„stp + 7r(A„)|5a„st|^)(ri, . . . , r„, Ai, . . . , A„) 

>e\Trn{uJ,Tn)+'K{An)) 

for some e > 0. Then, using the same techniques as in the proofs of Proposi- 
tions 3.2 and 3.1, one shows that the nondegeneracy (2.11) condition holds 
true. □ 

Remark 3.4. Note that the nondegeneracy condition holds true for 
every n (including n = 2) because we may use the noises coming from the 
jump times and the jump amplitudes at the same time. 



4. Numerical results. 
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4.1. Malliavin estimator. In this section, we compute the Deha of two 
European options: cah option with payoff <^(x) = (x — K)^ and digital option 
with payoff <j){x) = 1x>k- The asset {St)t>o follows a one-dimensional pure 
jump diffusion process. We use the notation from the beginning of Section 
3. 

We deal with two different pure jump diffusion models. The first is a 
Vasicek-type model, 

rt Jt 

(4.1) St = x- r{Su-a)du + y2aAi, 

•^0 i=l 

and the second is a geometrical model, 

(4.2) St = x+ / rSudu + ay2Sj,-Ai. 

^0 i=i ' 

In both models, we take Aj ~ AA(0, 1), i>l. That is, for alH > 1, Aj has the 
density p{x) = -^e^^^^ with p{x) = — \- Note that even if p is not bounded 
on M, the integration by parts formula holds by a truncation argument. 

Our aim is to compute dx^{^{^STy) using the integration by parts formula 
derived in the previous sections. We write 

9,.E(</.(5t))=E(0'(5t)5,5t) 

oo 

= E(0'(5t) a,5Tl{j^=o}) + E E(</''(5t) dxST\j^=n})- 

n=l 

For n > 1, we use the integration by parts formula on { = n} and obtain 

E((/)'(5t) 9,5Tl{J^=n}) = E(<^(5T)i^nl{J^=n}), 

where is a weight involving Malliavin derivatives of St and dxSx- Sum- 
ming over n = 1, 2, . . . , we obtain 

9,E(</>(5t)) =E(0'(5t)5.5t1{j^=o}) +E(</>(5T)i/j^(5T,9,5T)l{j^>i}). 

In order to compute the two terms in the right-hand side of the above 
equality, we proceed as follows. On {Jt = 0}, there is no jump on ]0, T], thus 
St and SxSt solve some deterministic integral equation. In the examples 
that we considered in this paper, the solutions of these equations are explicit, 
so this term is explicitly known. We may use the finite difference method. 
For the computation of the second term, we use a Monte Carlo algorithm. 
We simulate a sample ((r^)„gN> {^^)n&^), k = 1,. . . ,M, of the times and 
the amplitudes of the jumps and we compute the corresponding J^, 5^ and 
Hjf. . We then write 

, M 

E(<A(5T)ifj,(5T,a.5T)l{j,>i})=^^E<^('5r)^4l{4>i}- 
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We now compute the Malliavin weights Hj^{Sj^,dxS!p) for our examples. 
One may use Lemma 3.1, but in the particular cases that we discuss here, 
we have explicit solutions, so direct computations are much easier. 

• We first study the diffusion process defined by (4.1). We have the fol- 
lowing explicit expression for St on {Jt = n}: 

n 

(4.3) St = xe-^^ + a(l - e"^^) + aJ2 A.-e'^^^-^^). 

j=i 

We may use integration by parts with respect to the jump amplitudes, to 
the jump times or to both of them. 

* Jump amplitudes: Hj^^ has been calculated in [1]. Since Aj is Gaussian 
distributed for all i, the weight is tt{u>, Aj) = 1 and on {Jt = n}, 



(4.4) HniST,d,ST) 



* Jump times: suppose that n > 4 and Jt = ti. We use the weights 



Tri{uj,T.i) = {Ti+i - T,Y{Ti - Ti_iY and we have vr^ = a5l7^ S^-^di+i - 5i), 



where 6i = Ti — Ti-i. Then 

DiST = ciA.re-''^^-^') 

and 

LiST = -arAie-'^'^^-'^^\rTri + a{Si+i6ir-\6,+i-6i)), 

n 

1=1 

We define 

Bj = Afe^*-^^ [2rTTj + a{5j+i5jT-^{5j+i - 5j)]. 

Then 

DjCTsT = {arfe-'^''^{Aj_i5j-i - Aj+i5j+2 + Bj). 

Moreover, OxSt = e~^'^ so that Di OxSt = for all z = 1, . . . , n. 

We now have the expression for all of the terms involved in if„ and we 
obtain 

J^n — 

(4.5) 

ELi vrjAje^^«(^,_i5i_i - + Bi) 
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where 6- = E?=i 7r^Afe2^^^ 

For n = 1, 2, 3, we use integration by parts with respect to Ai only. Similar 
computations then give 



n — A • 

crAi 

• We study the jump diffusion process defined by (4.2). On { = n}, we 
have 

n 

5T = xe^^[](l + aAj). 

j=i 

We may not use integration by parts with respect to the jump times be- 
cause St depends on Ti, . . . , by means of Jt only. Therefore, we use the 
integration by parts formula with respect to the jump amplitudes only. On 
{Jt = n}^ the Malliavin weight is, in this case (see [1]), 

(4.6) Hn{ST,d,ST) = ^ + - 



where 



1 "A- 

^- = g(rT^' ''- = g(TT4, 



and 



-^(l + aA 



5. Numerical experiments. In this section, we present several numerical 
experiments in order to compare the Malliavin approach to the finite dif- 
ference method. We use the geometrical model and the Vasicek-type model 
and, in the latter, we use the Malliavin calculus with respect to the am- 
plitudes of the jump and to the jump times. We also look at two types of 
payoff: call options and digital options. 

The comparison is illustrated by some graphs (see Figures 1-3) on one 
hand and by empirical variance tables (see Tables 1-3) on the other. In 
Tables 1-3, we give the empirical variance of the two estimators denoted 
Var Mall and Var Diff and we compare them. We also include in our tables 
the value of the volatility a that we use and the corresponding variance of 
the underlying, denoted by Variance{St). We choose the parameter a in the 
following way. 



30 V. BALLY, M.-P. BAVOUZET AND M. MESSAOUD 

Delta of a Call European Oplion Eslimator , Vasicek model. K=S0=1 00,T=5.r=0,1,XHl.a=50 
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Fig. 1. Vasicek-type model. Delta of a European call option using MalUavm calculus 
based on jump times, on jump amplitudes and finite difference method. 



For the geometrical model, the variance of St is 



2 2rt/ a''\t 



Variance{St) = x e ^ {e' 



!)• 



Table 1 

Variance of the Malliavin JT estimator, AJ estimator and of the FD for call option in 

the Vasicek-type model 



cr 


Variance(ST) 


Var MallJT 


VarMallAJ 


Var Diff 


15.8114 


796.241 


0.0285123 


0.0106426 


0.0300379 


16.6667 


897.577 


0.0417219 


0.0115955 


0.0298567 


17.6777 


991.453 


0.0400695 


0.013123 


0.0298904 


18.8982 


1134.11 


0.0410136 


0.0144516 


0.0299574 


20.4124 


1313.42 


0.0433065 


0.0162378 


0.029862 


22.3607 


1584.9 


0.0400481 


0.0178726 


0.0298987 


25 


1967.53 


0.0407136 


0.0202055 


0.0299007 


28.8675 


2604.22 


0.0362728 


0.0224265 


0.0299651 


35.3553 


3961.31 


0.0343158 


0.0253757 


0.0297775 


50 


7890.4 


0.0333298 


0.0287716 


0.0299749 
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□etta of a Digital Europaan Option Estimator , Vasicek model. K=S0=1 M ,T=S.r=0.1 .0^5(1 
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Fig. 2. Vasicek-type model. Delta of an European digital option using Malliavm calculus 
based on the jump amplitudes, on the jump times and finite difference method. 



We take A = 1, r = 0.1, T = 5 and x = 100. Then for a G [0.1, 0.6], 1393.69 < 
Variance{ST) < 137264. 



Table 2 

Vasicek-type model. Variance of the Malliavin JT estimator, 
A J estimator and of the FD for digital option 



a 


Variancei^ST) 


Var MallJT 


VarMallAJ 


Var Diff 


15.8114 


796.241 


0.00144622 


7.18878e-5 


0.00514743 


16.6667 


897.577 


0.00254652 


7.3629e-5 


0.00459619 


17.6777 


991.453 


0.0018011 


7.85552e-5 


0.00496369 


18.8982 


1134.11 


0.0109864 


8.14005e-5 


0.00477995 


20.4124 


1313.42 


0.00177648 


8.1627e-5 


0.00386111 


22.3607 


1584.9 


0.00152777 


8.06193e-5 


0.00496369 


25 


1967.53 


0.0013786 


7.94341e-5 


0.0062497 


28.8675 


2604.22 


0.00100181 


7.5835e-5 


0.00551488 


35.3553 


3961.31 


0.000617271 


6.95225e-5 


0.00459619 


50 


7890.4 


0.000373802 


5.64325e-5 


0.00533116 
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Dfliha of a call European Option Estimator , Oeomalnc modal, K=S0=100,T=5.r=0.1,)l=|.(ji=n.2 
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Fig. 3. Geometrical model. Delta of an European digital option using MalUavm calculus 
based on the jump amplitudes and finite difference method. 



• For the Vasicek-type model, we have 

\ 2 

Variance(St) = 2aer'^''Hx - a) + — (1 - e'^^-*). 

2r 

We take A = 1, r = 0.1, T = 5, a = 10 and x = 100. Then for a G [16, 50] 
(note that the practitioners use o" = 20 to 30 in the short-term rate mod- 
ehng), 1471.3 < Variance^Sr) < 8563.69. 



Table 3 

Variance of the Malliavin estimator of the Delta and variance of the FD for digital option 



a 


Variancei^ST) 


Var Mall 


Var Diff 


0.1 


1405.06 


0.000263425 


0.00102718 


0.2 


6183.72 


0.000917207 


0.00164801 


0.3 


16005.5 


0.000885212 


0.00117345 


0.4 


42590.8 


0.000685313 


0.0013196 


0.5 


69018.7 


0.000531118 


0.000917399 


0.6 


130425 


0.000310461 


0.0003307 
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In all of our simulations, we have used a variance reduction method based 
on localization (analogous to the one introduced in [13] and [12]). We use 
the following abbreviations: 

• AJ: amplitude of the jumps. 

• JT: jump times. 

• FD: finite differences. 

• G: geometrical model. 

• V: Vasicek-type model. 

• Call: call option. 

• Dig: digital option. 

(G/Dig/AJ) then means that we are in the geometrical model (G) with a 
digital option (Dig) and we use an algorithm based on the amplitudes of the 
jumps (AJ). (G/Dig/AJ) versus (G/Dig/FD) means that we compare these 
two estimators. 

5.1. Numerical results for the Vasicek-type model. Here, we compare the 
Delta of European call and digital options obtained by using Malliavin cal- 
culus on the one hand and finite difference method on the other hand. 

• (V/Cah/AJ) versus (V/Call/JT) versus (V/Call/FD). 

• (V/Dig/AJ) versus (V/Dig/JT) versus (V/Dig/FD). 

In the call options case, both the graph and the variance table show that 
the Malliavin estimators and the finite difference variances are very close. 
The methods then give comparable results. 

In the digital options case, the algorithm based on Malliavin calculus is 
significantly better than the one based on finite difference. We also note that 
the Malliavin estimator using jump amplitudes gives less variance than the 
one based on jump times. 

5.2. Numerical results for the geometrical model. These are similar to 
those of the Vasicek-type model, as we can see for digital options. 

(G/Dig/AJ) versus (G/Dig/FD). 

5.3. Conclusions. 

• For a smooth payoff (as the call), the algorithms based on the Malliavin 
calculus (with respect to the jump times or amplitudes) give comparable 
results to those based on the finite difference method. 

• For a discontinuous payoff (as in the digital options), the algorithms based 
on Malliavin calculus give significantly better results than those based on 
the finite difference method. 
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